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In  this  paper,  the  impact  of  the  design  and  operating  parameters  of  a  tubular  Solid  Oxide  Fuel  Cell  (SOFC) 
is  studied  using  a  well-validated  steady-state  model.  The  profiles  of  species  concentrations  and  pressure 
in  the  flow  channels  when  the  cell  terminal  voltage  is  changed  are  studied.  In  addition,  both  the  axial  and 
radial  profiles  of  the  species  concentrations  inside  the  electrodes  are  presented.  The  model  is  also  used 
to  study  the  effects  of  the  parameters  that  can  significantly  influence  the  design  criteria  of  an  anode- 
supported  tubular  SOFC.  The  effects  of  the  flowrate  of  H2,  inlet  pressure,  and  the  cell  temperature  on 
the  power  output  from  the  cell  are  studied.  Cell  characteristic  parameters  such  as  the  porosity  of  the 
electrodes,  effective  diffusivity  of  the  species,  and  the  rate  of  the  electrochemical  reactions  are  varied  and 
their  impact  on  the  cell  performance  is  observed.  The  influence  of  the  cell  design  parameters  such  as  the 
thickness  of  the  electrodes  and  the  electrolyte  on  the  steady-state  polarization  curve  are  also  studied. 
Finally,  a  dimensional  study  is  presented.  In  the  dimensional  study,  the  radius  of  the  anode  flow  channel, 
the  length  of  the  cell,  and  the  annulus  size  are  varied  keeping  the  solid  volume  and  the  total  cell  volume 
constant.  This  study  shows  that  it  is  possible  to  design  guidelines  for  the  optimum  performance  of  the 
cell. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Efficient  conversion  of  the  fossil  fuels  to  electric  power  is 
absolutely  necessary  for  sustaining  the  growth  of  the  modern  civi¬ 
lization.  Solid  Oxide  Fuel  Cell  (SOFC)  is  a  high-temperature  fuel  cell 
that  is  envisaged  as  an  effective  mean  to  achieve  high  efficiency  in 
the  conversion  process.  For  commercialization  of  the  SOFC,  the  key 
focus  is  on  efficiency  and  the  cost.  This  paper  presents  the  effects  of 
operating  conditions  of  the  cell,  its  characteristic  parameters,  and 
the  design  parameters  on  the  cell  performance.  Such  a  study  can 
be  extremely  useful  in  deciding  the  decision  variables  for  design¬ 
ing  an  optimum  cell  and  for  determining  the  optimum  operating 
conditions  of  a  cell. 

Studies  on  transport  fields,  especially  concentration  and  pres¬ 
sure  fields,  can  reveal  important  information  about  the  cell  such  as 
gas  starvation  that  needs  to  be  considered  during  the  cell  design. 
However  such  studies  are  minimal  in  the  existing  literature  [  1  -3 ].  A 
few  studies  can  be  found  that  have  discussed  the  transport  profiles 
only  in  the  main  flow  direction  [4-6].  However,  significant  gradients 
may  develop  in  other  directions  too.  For  example,  the  radial  gradi- 
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ent  of  concentration  inside  the  electrodes  of  a  tubular  SOFC  can  be 
significant  especially  at  high  current  density.  Therefore  important 
observations  can  be  made  by  analyzing  the  transport  profiles  in  the 
radial  direction  of  a  tubular  SOFC.  This  analysis  can  also  indicate 
when  it  is  important  to  consider  these  phenomena  for  a  predic¬ 
tive  model.  Several  authors  [7-9]  have  studied  the  reactant  and 
product  concentration  profiles  in  the  directions  other  than  the  bulk 
flow  direction.  However,  these  results  have  been  developed  for  pla¬ 
nar  SOFCs.  Also,  the  concentration  profiles  have  been  reported  for 
a  particular  operating  condition.  Studies  on  the  profiles  of  the  dif¬ 
ferent  transport  fields  of  the  cell  at  different  operating  conditions 
are  rare. 

Significant  changes  in  the  cell  performance  can  result  if  the 
flow  rates  of  the  reactants  are  changed.  For  a  pressure-driven  flow, 
change  in  flow  causes  change  in  the  pressure  in  the  flow  channels. 
However,  pressure  changes  may  not  be  significant  for  gaseous  reac¬ 
tants  at  higher  temperatures.  On  the  other  hand,  the  cell  operating 
pressure  can  be  changed  at  the  same  flowrate  by  installing  pres¬ 
sure  control  valves.  In  this  work,  the  influence  of  H2  flowrate  and 
the  operating  pressure  on  the  cell  performance  are  studied. 

The  impact  of  the  oxidant  flowrate  on  the  cell  performance  is 
reported  in  the  work  of  Jia  et  al.  [5].  In  their  work,  the  cell  per¬ 
formance  decreases  slightly  as  02  flow  is  increased.  The  reason  is 
the  decrease  in  the  cell  temperature  due  to  the  flow  of  excess  air. 
However,  the  change  in  the  temperature  depends  upon  the  control 
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Nomenclature 

Aj PB  active  area  for  electrochemical  reactions  (m2) 

C  concentration  of  species  (mol  m-3 ) 

Dj  eff  effective  diffusivity  of  species  j  (m2  s-1 ) 

Dj_n  binary  diffusion  coefficient  between  species  j  and  n 

(m2  s-1) 

Dj  K,e ff  effective  Knudsen  diffusivity  for  species  j  (m2  s-1 ) 

Ii  current  generated  by  the  ith  CV  (A) 

rin, ac  inside  radius  of  the  anode  channel  (m) 
rin,cc  inside  radius  of  the  cathode  channel  (m) 
rout ,cc  outside  radius  of  the  cathode  channel  (m) 

uac  axial  velocity  in  the  anode  channel  (m  s-1 ) 
vac  radial  velocity  in  the  cathode  channel  (m  s-1 ) 

Greek  symbols 
8  porosity 

Subscripts 

ac  anode  channel 

an  anode 

cc  cathode  channel 

ca  cathode 

i  index  for  the  control  volume 

j  index  for  the  species:  H2,  02,  N2,  H20 


strategy  adopted  for  maintaining  the  temperature  of  the  cell.  For 
the  cell  being  studied,  as  the  cell  temperature  is  being  maintained 
by  a  feedback  controller,  an  increase  of  the  heat  loss  because  of 
the  excess  air  will  be  compensated  by  the  heating  coils.  The  air¬ 
flow  is  maintained  constant  for  this  cell.  The  effects  of  changes 
in  the  H2  flowrate  are  observed  in  this  work.  Such  a  study  can  be 
found  in  the  work  of  Costamagna  and  Honegger  [10].  However  the 
effects  of  increasing  H2  flowrate  at  different  cell  temperatures  are 
not  reported. 

A  number  of  authors  have  discussed  the  influence  of  pressure 
[1,4,3].  It  is  observed  in  the  study  of  Ni  et  al.  [3]  that  the  current 
density  of  an  anode-supported  cell  increases  from  1.5  A  cm-2  to 
6  A  cm-2  as  the  cell  pressure  increases  from  0.5  bar  to  5  bar.  In  their 
work,  the  transport  of  the  gaseous  species  inside  the  flow  chan¬ 
nels  is  not  considered  at  all.  Also,  a  simplified  model  is  considered 
for  the  concentration  overpotential  in  the  electrodes.  Such  models 
may  not  truly  capture  the  rise  in  the  concentration  overpoten¬ 
tial  as  the  current  density  increases.  For  a  planar  anode-supported 
SOFC,  the  work  of  Hussain  et  al.  [1]  shows  that  the  %  increase  in 
the  cell  performance  is  gradually  reduced  as  the  cell  pressure  is 
increased.  In  their  work,  the  observed  current  density  at  0.5  V  is 
about  0.8  A  cm-2, 1.05  A  cm-2,  and  1.2  A  cm-2  fori  atm,  3  atm,  and 
5  atm,  respectively.  In  our  work,  similar  studies  are  done  for  a  tubu¬ 
lar  anode-supported  SOFC  by  increasing  the  inlet  pressures  of  the 
gas  flow  channels. 

The  influence  of  temperature  on  the  cell  performance  can  be 
seen  in  the  work  of  Hussain  et  al.  [1];  Jia  et  al.  [5];  Nagata  et  al. 
[11];  Ni  et  al.  [3].  Ni  et  al.  [3]  have  shown  the  performance  to  be 
0.2  A  cm-2, 1 A  cm-2,  and  3  A  cm-2,  respectively,  at  600  °C,  800  °C, 
and  1000  °C,  respectively  at  0.5  V.  A  significant  gain  in  the  cur¬ 
rent  density  at  a  given  voltage  is  expected  as  the  temperature  is 
increased  mainly  because  of  the  increase  in  the  conductivity  of  the 
electrolyte  and  the  reduction  of  the  activation  overpotential.  How¬ 
ever  the  magnitude  of  gain  reported  in  the  work  of  Ni  et  al.  [3]  is 
unlikely  and  may  be  due  to  a  simplified  model  for  the  concentration 
overpotential  as  mentioned  before.  The  simulation  results  of  Hus¬ 
sain  et  al.  [1  ]  show  the  current  density  to  change  from  0.75  A  cm-2 
at  800  °C  to  1 .2  A  cm-2  at  1000  °C.  However  these  studies  are  done 


for  a  planar  SOFC.  The  work  of  Jia  et  al.  [  5  ]  is  on  a  cathode-supported 
tubular  SOFC.  An  insignificant  gain  in  the  cell  performance  with  an 
increase  in  the  temperature  is  observed  in  their  work.  The  reason 
may  be  a  comparably  low  increase  in  the  temperature  of  PEN  due  to 
increase  of  the  reactant  inlet  temperature.  It  can  be  seen  from  the 
above  literature  review  that  the  gain  predicted  by  a  model  due  to 
the  increase  in  the  temperature  varies  widely.  It  not  only  depends 
on  the  phenomena  considered  in  a  model,  but  also  can  vary  widely 
based  on  the  cell  configuration  and  the  control  strategy  adopted. 
In  this  work,  studies  on  the  effects  of  the  cell  temperature  on  an 
anode-supported  tubular  SOFC  are  presented. 

Porosity  can  have  a  strong  effect  on  the  concentration  overpo¬ 
tential,  activation  overpotential,  and  the  ohmic  losses  of  the  cell. 
Significant  differences  in  observations  can  be  seen  in  different 
studies  [1,3].  Hussain  et  al.  [1]  report  the  effect  of  changes  in  the 
porosity  in  the  range  of  0.3-0.7.  In  this  range,  the  cell  performance 
drops  monotonically  as  the  porosity  increases.  Ni  et  al.  [3]  iden¬ 
tify  an  optimal  porosity  of  0.4  by  studying  different  porosities  in 
a  range  of  0.2-0.6.  In  the  present  work,  the  effects  of  the  porosi¬ 
ties  of  the  electrodes  are  considered  separately.  Optimum  porosity 
for  both  the  electrodes  are  also  identified  for  a  given  operating 
condition. 

Effective  diffusivity  of  the  reactants  can  change  because  of  sev¬ 
eral  reasons,  such  as  porosity,  tortuousity,  temperature,  and  pore 
radius.  Although  all  of  them  affect  effective  diffusivity,  their  effect 
on  other  loss  mechanisms  may  be  different.  Instead  of  studying 
any  specific  reason  for  the  increase  in  effective  diffusivity,  a  study 
is  done  by  varying  the  effective  diffusivity  in  general.  Effective 
diffusivities  of  the  species  in  the  anode  and  cathode  are  varied 
separately.  Significant  differences  are  observed  when  the  effective 
diffusivity  is  changed  for  the  species  in  the  cathode  compared  to 
that  of  the  anode. 

The  rate  constants  of  the  electrochemical  reactions  affect  the 
activation  overpotential  of  the  cell  strongly.  Rate  constants  can 
change  based  on  the  microstructural  properties  such  as  the  grain 
size,  contact  angle  between  electron  and  ion  conducting  particles, 
and  grain  boundary  effects  [12,13].  Instead  of  changing  these  fac¬ 
tors  separately,  the  effects  of  changes  in  the  overall  rate  constant 
on  the  performance  of  the  cell  are  observed.  Both  the  anodic  and 
cathodic  rate  constants  are  changed  and  significant  differences  in 
their  effects  are  observed.  Studies  on  the  effects  of  the  effective 
diffusivities  and  rate  constants  are  minimal  in  the  literature. 

The  cell  design  parameters  can  have  a  strong  effect  on  the  cell 
performance.  The  thickness  of  the  electrodes  can  significantly  influ¬ 
ence  both  the  concentration  overpotential  and  the  ohmic  losses  of 
a  cell.  The  thickness  of  the  anode  of  an  anode-supported  tubular 
SOFC  will  affect  the  TPB  surface  area  for  both  the  anodic  and  the 
cathodic  reactions  for  a  given  radius  of  the  anode  channel.  In  con¬ 
trast,  the  cathode  thickness  does  not  have  any  effect  on  the  TPB 
surface  area  of  both  the  anode  and  the  cathode  for  a  LSM  cathode 
as  the  active  thickness  is  not  expected  to  be  beyond  10-100  p,m  [14]. 
Interesting  differences  in  the  cell  performance  are  observed  as  the 
thicknesses  of  the  electrodes  are  changed.  Increasing  the  thickness 
of  the  electrolyte  can  change  the  ohmic  resistance  in  the  electrolyte 
significantly  and  bring  down  the  performance  to  a  great  extent. 
However  the  change  in  the  thickness  of  the  electrolyte  also  affects 
the  current  path  length  in  the  cathode  and  the  active  surface  area 
for  the  cathodic  reaction. 

Although  a  number  of  studies  exists  on  the  system  level  opti¬ 
mization  of  SOFC  [15-18],  studies  on  unit  level  optimization  of 
SOFC  are  rare.  One  of  the  main  challenges  for  the  commercializa¬ 
tion  of  SOFC  is  the  material  cost.  At  the  same  time,  volume  packing 
of  a  SOFC  is  very  important  for  both  man-portable  and  stationary 
powerhouse  applications.  To  optimize  the  amount  of  solid  material 
and  the  volume  packing,  an  optimal  combination  of  length,  anode 
channel  radius,  and  annulus  size  is  required  for  a  SOFC.  In  order  to 
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identify  whether  such  an  optimal  combination  exists,  a  dimensional 
study  is  presented. 

The  key  aspects  of  this  paper  are  as  follows: 

(1)  To  get  representative  results  from  the  sensitivity  study  of  a 
cell,  the  used  model  should  be  validated  with  experimental 
results  over  a  broad  operating  range  as  is  expected  during  the 
sensitivity  study.  If  the  model  itself  is  not  valid  in  a  certain 
domain,  the  sensitivity  results  from  such  a  model  will  have  lit¬ 
tle  significance.  The  studies  presented  in  this  paper  have  been 
done  by  a  well-validated  model.  During  the  validation  study 
[19],  we  have  identified  the  domain  of  the  operating  condi¬ 
tions  in  which  the  model  is  valid.  The  sensitivity  studies  have 
been  done  in  this  valid  domain.  Some  of  the  sensitivity  stud¬ 
ies  presented  in  this  work  have  been  done  for  the  planar  cell 
for  various  authors,  but  are  rare  for  anode-supported  tubu¬ 
lar  SOFC.  Additional  sensitivity  and  dimensional  studies  are 
done  to  provide  insight  into  the  design  of  the  cell,  its  operat¬ 
ing  condition,  and  improvement  in  the  cell  efficiency  by  doing 
dimensional  optimization.  The  main  focus  in  all  the  studies  has 
been  on  the  aspects  that  can  improve  the  cell  performance. 

(2)  This  study  shows  the  profiles  of  concentration  and  pressure 
both  in  the  radial  and  the  axial  directions  inside  the  flow  chan¬ 
nels  and  the  electrodes  of  an  anode-supported  tubular  SOFC. 
The  study  presents  the  limiting  phenomena  at  various  operat¬ 
ing  conditions  and  possibility  of  gas  starvation  at  high  current 
density.  This  suggests  that  consideration  of  the  mass  transfer 
is  very  important  for  a  predictive  model  over  a  broad  operating 
range. 

(3)  Effects  of  change  in  the  flowrate  of  H2  at  different  cell  temper¬ 
atures  are  also  studied.  Such  a  study  is  important  for  deciding 
the  optimum  flowrate  of  H2  as  the  cell  temperature  is  varied. 

(4)  The  effects  of  the  operating  pressure  on  this  tubular  cell  have 
been  presented.  The  study  suggests  that  the  operating  pressure 
should  be  determined  based  on  the  operating  condition  of  the 
cell  and  on  its  hardware  design.  Else,  the  parasitic  loss  due  to 
compression  can  decrease  the  overall  system  efficiency. 

(5)  A  study  is  done  to  show  the  effects  of  the  temperature  on 
the  cell  performance.  For  this  anode-supported  cell  of  thin- 
electrolyte  configuration,  it  is  observed  that  the  improvement 
in  the  cell  performance  is  minimal,  if  at  all,  beyond  a  cer¬ 
tain  temperature.  Because  of  the  parasitic  loss  for  maintaining 
higher  temperature,  the  overall  system  efficiency  may  actu¬ 
ally  decrease  beyond  certain  temperature.  The  study  can  be 
useful  for  deciding  the  optimum  operating  temperature  for  a 
particular  cell. 

(6)  The  effects  of  the  porosity  of  the  electrodes  on  the  cell  perfor¬ 
mance  are  also  studied.  Significant  differences  in  the  effects  are 
observed  due  to  change  in  the  porosity  of  the  anode  and  the 
cathode  for  similar  changes.  The  work  shows  the  necessity  of 
an  optimum  porosity  for  getting  improved  performance  from 
a  cell. 

(7)  Because  of  the  near-atmospheric  operating  pressure  at  the 
cathode  channel  of  the  cell  and  due  to  lower  diffusivity  of  02 
than  H2,  the  change  in  the  effective  diffusivity  of  02  in  the 
cathode  is  found  to  have  significant  effect  in  the  concentration 
polarization  of  the  cell  compared  to  that  of  H2  in  the  anode. 
These  aspects  need  to  be  considered  at  the  design  stage  itself 
for  an  optimized  cell. 

(8)  The  study  on  the  effects  of  the  reaction  rate  is  useful  for 
designing  the  electro-catalysts.  The  study  shows  the  operat¬ 
ing  regions  where  considerable  improvement  is  possible  by 
enhancing  the  rate  of  reactions. 

(9)  This  work  also  presents  a  study  by  changing  the  thickness 
of  the  electrodes  and  the  electrolytes.  A  change  in  the  thick¬ 
ness  of  a  particular  cell  component  can  be  decided  based  on 
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the  design  objective.  For  example,  a  decrease  in  the  anode 
thickness  improves  the  cell  performance  still  the  cell  hits  the 
limiting  current  density  which  gets  decreased  due  to  decrease 
in  the  anode  thickness.  So  the  best  anode  thickness  can  be  cho¬ 
sen  based  on  a  design  philosophy  satisfying  the  mechanical 
integrity  of  the  cell  and  constraints  of  manufacturing. 

(10)  The  dimensions  of  a  SOFC,  particularly  that  of  a  cylindrical  cell, 
has  a  significant  effect  on  the  performance  of  the  cell.  If  the 
volume  of  the  PEN  is  kept  constant,  a  decrease  in  the  radius 
of  the  anode  channel  will  result  in  a  longer  length  of  the  cell. 
Again  if  the  total  volume  of  the  cell  is  constant  for  an  anode- 
supported  tubular  SOFC,  the  change  in  the  radius  of  the  anode 
channel  will  result  in  a  different  size  of  the  annulus  of  the 
cathode  channel.  As  the  cell  length,  the  radius  of  the  anode 
channel,  and  the  annulus  size  of  the  cathode  channel  change, 
the  pressure  drop  across  the  flow  channels  changes  along  with 
other  transport  fields.  For  a  pressure-driven  flow,  the  maxi¬ 
mum  pressure  is  expected  at  the  inlet  of  the  gas  channel.  So,  if 
the  cell  is  run  at  the  maximum  pressure  that  the  cell  can  with¬ 
stand,  then  an  increased  pressure  drop  in  the  flow  channels 
will  cause  an  increase  in  the  activation  and  concentration  over¬ 
potential  of  the  cell.  Also  the  electrons  have  to  flow  from/to 
the  current  collectors  to/from  the  triple  phase  boundary  (TPB, 
where  the  electrochemical  reactions  take  place)  because  of  the 
electrochemical  reactions.  As  the  current  collectors  are  typ¬ 
ically  placed  only  over  a  small  region  of  the  circumference, 
the  current  path  length  changes  as  the  cell  radius  changes 
affecting  the  ohmic  resistance  of  the  cell.  The  change  in  the 
annulus  size  of  the  cathode  channel  affects  the  cathode  acti¬ 
vation  and  concentration  polarization  which  is  significant  for  a 
tubular  SOFC  especially  at  a  high  current  density.  So  a  decrease 
in  the  radius  of  the  anode  channel  for  an  anode-supported 
tubular  SOFC  will  increase  pressure  drop  in  the  anode  chan¬ 
nel  because  of  the  increased  cell  length  and  will  decrease  the 
ohmic  resistance  because  of  a  shorter  current  path  (assum¬ 
ing  the  thickness  of  the  electrodes  and  electrolytes  remains 
unchanged).  These  competing  loss  mechanisms  suggest  the 
existence  of  an  optimum  combination  of  the  dimensions  of 
the  anode  channel  radius,  cell  length,  and  the  cathode  chan¬ 
nel  annulus.  A  dimensional  optimization  study  shows  that  the 
optimum  power  density  is  achievable  at  an  inner  radius  of  the 
anode  that  can  be  manufactured  commercially. 

2.  Background  on  the  experimental  SOFC  and  the 
optimization  oriented  model 

The  schematic  representation  of  the  anode-supported  tubular 
industrial  cell  used  in  this  study  is  shown  in  Fig.  1.  The  cell  is 
of  counterflow  configuration.  As  H2  flows  down  the  anode  chan- 
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Fig.  2.  Illustration  of  the  coordinate  system. 


nel,  it  moves  to  the  anode  TPB  for  the  electrochemical  reaction. 
H20,  product  of  the  reaction,  flows  from  the  TPB  to  the  anode  flow 
channel.  In  the  cathode  channel,  air  flows  in  the  opposite  direction 
to  that  of  H2  in  the  anode  channel.  02  flows  to  the  cathode  TPB  for 
taking  part  in  the  electrochemical  reaction.  A  detailed  steady-state 
mode  for  this  cell  has  been  presented  in  our  previous  works  [19,20]. 
The  model  was  also  extensively  validated  in  those  publications  at 
several  different  reactant  flowrates  and  temperatures.  In  what  fol¬ 
lows,  we  briefly  describe  the  major  assumptions  in  the  model,  the 
governing  equations,  the  key  features  of  the  steady-state  model,  and 
the  boundary  conditions.  The  following  description  of  the  model  is 
provided  as  a  quick  reference  for  the  study  presented  later  in  this 
work.  For  all  the  modeling  equations  and  detailed  validation  results, 
an  interested  reader  is  referred  to  our  previous  works  [19,20]. 

2  A.  Model  assumptions 


significant  if  the  current  density  is  high.  As  the  intended  use  of 
this  model  is  to  maximize  the  performance  of  the  cell  by  opti¬ 
mizing  the  cell  dimensions,  it  is  required  to  consider  both  the 
axial  and  radial  variations  in  the  concentration  fields.  The  species 
conservation  equation  in  the  anode  channel  is  written  as  [19]: 

dCj,ac  1  d(Cj,acrVac)  d{Cj,acUac)  _  1  3  (  dCj  ac\  _  „ 

3 t  r  dr  3 z  ij  r  dr  y  dr  ) 

(1) 


where  j  =  H2 ,  H2  0,  vac  is  the  r-component  of  velocity  in  the  anode 
channel  and  Dy  is  the  binary  diffusivity  of  H2-H20  system.  Similar 
equations  are  considered  for  the  components  02  and  N2  in  the 
cathode  channel. 

•  Both  the  r-component  and  z-component  of  the  equation  of 
motion  are  considered  in  the  gas  flow  channels.  z-Component 
of  the  equation  of  motion  in  the  anode  channel  is  written  as  [  19] : 
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dt 


■  +  Va< 


3 Uac  dUac  \ 

lk+Uaclk  ) 


dPac  Pac  3  /  3 Uac  \ 

rdr  \  ~dF  J 
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where  /xac  is  the  viscosity  in  the  anode  channel.  Similarly,  the 
r-component  of  the  equation  of  motion  in  the  anode  channel  is 
written  as 
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Similar  equations  are  considered  in  the  cathode  channel  for  the 
conservation  of  momentum. 

•  Species  conservation  equations  are  considered  in  the  electrodes. 
In  the  anode,  the  species  conservation  equation  is  written  as 


The  following  assumptions  are  made  as  follows: 

•  The  flow  is  newtonian.  The  mixture  viscosity  is  calculated  assum¬ 
ing  that  two  polar  molecules  interact  according  to  Stockmayer 
potential. 

•  The  present  model  is  isothermal.  The  cell  is  kept  inside  an  well- 
insulated  furnace.  The  inner  wall  of  the  insulator  has  heating 
coils;  current  through  which  is  controlled  to  maintain  the  cell 
temperature  by  a  temperature  controller.  Due  to  this  config¬ 
uration  and  the  controller  setup,  the  spatial  variation  of  the 
temperature  along  the  cell  is  neglected.  Experimentally,  this  is 
found  to  be  generally  valid  for  a  single  industrial  cell  operating 
with  pure  H2. 

•  The  primary  mechanism  of  mass  transport  in  the  porous  elec¬ 
trodes  is  considered  to  be  diffusive.  In  other  words,  convective 
transport  of  mass  inside  porous  electrodes  is  considered  negligi¬ 
ble  compared  to  the  diffusive  transport  [21  ]. 

•  Control  volume  approach  is  adopted  [20,22]  in  this  model.  In  the 
control  volumes  of  the  gas  channels,  densities  ( pac  and  pcc)  are 
assumed  to  be  constant.  This  assumption  simplifies  the  momen¬ 
tum  balance  equations. 


2.2.  Governing  equations  and  key  features  of  the  steady-state 
model 

The  coordinate  system  and  the  nomenclature  followed  in  the 
radial  and  axial  directions  are  shown  in  Fig.  2. 

The  following  phenomena  are  included  in  the  model: 

•  Both  the  axial  and  radial  variations  of  the  species  concentration 
fields  are  considered  in  the  flow  channels.  It  is  shown  later  in  this 
study  that  the  radial  variations  in  the  concentration  field  may  be 


&an 


dCj,an 

dt 


-  Dj,e ff 


d2Cj,an\ 

3z2  J 


(4) 


where  san  is  the  porosity  of  the  anode,  Cjan  is  the  concentration 
of  species  j  in  the  anode,  Dj  eff  is  the  effective  diffusivity  of  species 
j  in  the  anode  and  j  =  H2,  FI20.  In  the  cathode,  the  conservation 
equations  are  considered  for  the  species  02  and  N2.  Both  binary 
and  Knudsen  diffusivity  are  considered.  The  calculation  of  the 
effective  diffusivity  can  be  found  in  Bhattacharyya  et  al.  [20]. 

•  Several  loss  mechanisms  in  the  cell  reduces  the  cell  terminal 
voltage  from  the  theoretical  Nernst  potential.  The  major  losses 
are  concentration,  activation,  and  ohmic  losses.  Concentration 
losses  take  place  because  of  the  flow  of  the  reactants/products 
through  the  electrodes  to/from  the  TPB.  The  concentration  losses 
are  accounted  for  by  considering  the  species  conservation  equa¬ 
tions  in  the  electrodes  mentioned  previously.  The  electrochemical 
reactions  are  considered  to  follow  Butler-Volmer  kinetics.  A  ‘sin/T 
approximation  of  the  Butler-Volmer  equation  is  done  for  captur¬ 
ing  the  activation  losses.  An  Arrhenius-type  expression  is  written 
for  the  exchange  current  density  [19].  The  model  presents  an 
equivalent  circuit  approach  for  the  calculation  of  the  ohmic  resis¬ 
tance  considering  the  current  path  length  in  the  electrodes. 

•  Detailed  models  employing  mass  and  momentum  conservation 
equations  along  with  the  electrochemical  reactions  are  typically 
solved  in  the  commercial  CFD  softwares  [23,24].  In  this  model,  a 
MAPLE-MATLAB  combination  is  employed  to  solve  the  system. 
Flexibility,  computational  issues,  and  limitations  of  this  approach 
have  been  discussed  elsewhere  [19,20].  As  the  intended  use  of  the 
model  is  to  do  optimization  studies  of  the  cell,  the  model  is  solved 
in  MATLAB  platform  so  that  the  powerful  optimization  algorithms 
available  in  MATLAB  can  be  utilized. 

•  The  model  has  been  validated  over  a  wide  operating  range  of 
operating  conditions.  In  the  process  of  validation,  important  phe- 
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Table  1 

Boundary  conditions. 


Number 

Location 

Boundary  condition 

Comments 

Anode  gas  flow  channel 

BC1  z  —  0,  Vr 

Cj,ac  =  Q,ac, inlet 

j  —  H2,  H20,  inlet  conditions 

BC2 

II 

o 

Jf 

3-¥=° 

No  flux 

BC3 

r  =  rin,ac>  VZ 

Flux  continuity 

BC4 

> 

o' 

II 

N 

Wac  =  Wac.inlet 

Inlet  condition 

BC5 

II 

O 

duac  _  n 
dr  ~U 

Force  balance 

BC6 

II 

Jf 

Uac  =  0 

No  slip 

Anode 

BC7 

r  =  Hn.ac,  VZ 

Cj,an  =  Cj,ac 

Continuity  at  surface,  j  =  H2,  H20 

BC8 

s- 

> 

o' 

II 

N 

-Dj,eff^=0 

Zero  flux  condition 

BC9 

z  =  L,  Vr 

1 

ii 

O 

Zero  flux  condition 

BC10 

r  =  routtan,  Vz 

=  £ 

Flux  to  TPB  equals  rate  of  consumption 

Cathode  gas  flow  channel 

BC11  z  =  L,  Vr 

E/,cc  =  Q,cc,  inlet 

Inlet  conditions,  j  =  02,  N2 

BC12 

r  =  rout,cc  Vr 

^££-0 
dr  ~  U 

Impermeability 

BC13 

r  =  rin,cc>  VZ 

-DV~W^  +vccCo2,cc  =  -Do2,e«3-^1 

Flux  continuity 

BC14 

II 

~Dij  +  ^ccCn2,cc  =  0 

Flux  continuity 

BC15 

z  —  L,  Vr 

Wcc  =  ^cc.inlet 

Inlet  condition 

BC16 

r  =  rinxc,  Vz 

Ucc  =0 

No  slip 

BC17 

r  =  rout, cc,  Vz 

Ucc  =0 

No  slip 

BC18 

z  =  L,  Vr 

VCc  =  0 

Inlet  condition 

BC19 

r  =  rout,  co  Vz 

l^cc  =0 

Impermeability 

BC20 

r  =  rin,co  VZ 

p.  3  C0  cc 

—Dij  — 97 —  —  vcc  Cn2,cc 

Flux  condition 

Cathode 

BC21 

r  =  rin,co  VZ 

P 

s 

II 

p 

Continuity  at  surface 

BC22 

z  =  0,  Vr 

n  3c0  ca 

,eff  97  =  0 

Zero  flux  condition 

BC23 

> 

-J 

II 

N 

n  3C0  ca 

-L>o2,eff— a*—  =  0 

Zero  flux  condition 

BC24 

r  =  rin,Ca,  VZ 

n  a  dc02,ca  It 

-Oo2,eff/iTPB  ,ca,i—dT~  = 

Flux  to  TPB  equals  rate  of  consumption 

nomena  in  different  operating  regions  have  been  identified  and 
their  effects  on  the  cell  performance  have  been  highlighted.  A 
detailed  account  of  the  validation  studies  can  be  found  in  Bhat¬ 
tacharyya  [19]. 

2.3.  Boundary  conditions 

Inlet  conditions,  continuity  at  interface,  impermeability,  no  slip 
and  continuity  of  flux  are  used  as  the  boundary  conditions  for  solv¬ 
ing  the  system  of  equations.  The  boundary  conditions  are  listed  in 
Table  1. 


3.  Studies  on  the  transport  profiles 

Studies  on  transport  profiles  are  done  for  the  cell  temperature 
of  850  °C  and  36  ml  min-1  H2  flowrate.  The  study  is  done  for  differ¬ 
ent  cell  terminal  voltages— 0.97  V,  0.85  V,  and  0.65  V.  Even  though 
significantly  different  profiles  are  expected  at  different  operating 
conditions,  this  particular  operating  condition  and  the  cell  termi¬ 
nal  voltages  are  chosen  as  interesting  system  characteristics  are 
observed  in  these  conditions. 


CO 
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Fig.  3.  Profile  of  H2  concentration  in  anode  channel. 
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Fig.  3  shows  the  profile  of  H2  concentration  along  the  flow 
channel.  The  decrease  of  the  concentration  in  the  axial  direction 
is  significantly  higher  than  that  in  the  radial  direction.  As  the  cell 
voltage  is  decreased,  more  H2  is  consumed  causing  a  steeper  fall  in 
the  concentration.  The  radial  variation  of  H2  concentration  is  not 
clear  in  Fig.  3  because  of  its  lower  magnitude  compared  to  that  of 
the  axial  variations  in  FI2  concentration.  Flowever,  as  the  voltage  is 
lowered,  the  radial  variation  can  be  significant.  To  make  it  clear,  the 
radial  variation  of  FI2  concentration  with  the  change  in  voltage  at  a 
particular  axial  location  is  shown  in  Fig.  4.  It  is  observed  that  FI2  con¬ 
centration  decreases  from  the  center  of  the  channel  (r  =  0)  towards 
the  wall.  The  radial  decrease  of  the  concentration  is  negligible  when 
the  cell  voltage  is  0.97  V.  However,  this  becomes  increasingly  sig¬ 
nificant  as  the  cell  voltage  is  lowered.  It  is  seen  in  both  Figs.  3  and  4 
that  the  concentration  of  H2  decreases  at  a  particular  location  along 
the  center  of  the  channel  with  a  decrease  in  cell  terminal  volt¬ 
age.  Because  of  the  cumulative  effects  of  both  the  radial  and  axial 
concentration  changes,  the  concentration  of  H2  in  anode/flow  chan¬ 
nel  interface  decreases  significantly  as  the  cell  terminal  voltage  is 
decreased.  These  findings  show  that  it  is  important  to  consider 
radial  variation  in  the  concentrations  fields  if  higher  currents  are 
generated  or  if  the  radius  is  used  as  an  optimization  variable. 


0.65  V 


0.85  V 


0.97  V 


02  I7non-dirne'1s'0^!"Ze^ 
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Fig.  5.  Profile  of  H20  concentration  in  anode  channel. 


Fig.  6.  Pressure  profile  in  axial  direction  in  anode  channel. 


Fig.  5  shows  the  concentration  profile  of  H20  in  the  anode  chan¬ 
nel.  As  H20  is  being  generated  due  to  the  electrochemical  reaction, 
concentration  keeps  increasing  along  the  length  of  the  anode  chan¬ 
nel.  Although  the  gain  in  H20  concentration  is  negligible  at  the 
cell  terminal  voltage  of  0.97  V,  it  increases  considerably  along  the 
length  of  the  anode  channel  as  the  cell  voltage  is  lowered.  As  before, 
the  change  in  the  concentration  of  H20  in  the  radial  direction  is  not 
clear  in  Fig.  5  as  it  is  insignificant  compared  to  the  axial  change 
in  the  concentration.  However,  a  close  observation  shows  a  radial 
increase  in  the  concentration  from  the  center  of  the  channel  to 
the  wall  and  this  increase  becomes  significant  as  the  cell  terminal 
voltage  is  lowered. 

At  this  operating  condition,  the  pressure  drop  in  both  radial 
and  axial  directions  inside  the  anode  channel  are  negligible.  Fig.  6 
shows  the  pressure  drop  in  the  anode  channel  in  the  axial  direction. 
Although  the  pressure  drop  is  negligible,  it  keeps  increasing  as  the 
voltage  decreases.  The  main  reason  for  the  increased  pressure  drop 
with  voltage  is  the  generation  of  more  water  with  the  decrease  in 
voltage.  Water  is  about  four  times  more  viscous  than  H2  at  the  same 
temperature.  Therefore  increased  water  vapor  concentration  down 
the  channel  increases  the  pressure  drop  as  the  cell  voltage  falls. 

Fig.  7  shows  the  concentration  profile  of  H2  in  the  anode.  In 
the  plot,  ‘radial  distance  =  0’  indicates  the  anode-flow  channel 
interface  and  ‘radial  distance  =  V  indicates  the  anode-electrolyte 
interface.  Again  it  is  observed  that  the  overall  fall  in  the  concentra¬ 
tion  of  H2  in  the  radial  direction  is  much  less  compared  to  that  in 
the  axial  direction  for  all  the  three  voltages  studied.  However,  the 
fall  in  the  concentration  is  found  to  increase  in  the  radial  direction 
as  the  voltage  is  decreased. 

In  the  cathode  flow  channel,  a  significant  drop  in  02  concentra¬ 
tion  in  the  radial  direction  is  observed  in  Fig.  8  as  the  cell  terminal 
voltage  is  decreased.  In  this  figure,  ‘radial  distance  =  0’  denotes  the 
cathode-flow  channel  interface  and  ‘radial  distance  =  V  denotes 
the  insulator-flow  channel  interface.  It  might  be  mentioned  that 
in  the  cathode  channel,  the  pressure  is  lower  than  that  of  the 
anode  channel  and  also  air  (unlike  pure  H2  in  anode  channel)  flows 
through  it.  At  0.65  V,  towards  the  end  of  the  cathode  channel,  the  02 
concentration  is  found  to  be  approaching  zero  at  the  cathode-flow 
channel  interface.  This  observation  strongly  suggests  the  necessity 
of  considering  the  radial  variation  of  02  concentration  inside  the 
cathode  channel  especially  for  a  cell  operating  at  near-atmospheric 
pressure  and  high  current  density. 
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Fig.  7.  Concentration  profile  of  H2  in  anode. 


The  pressure  drop  in  the  cathode  flow  channel  is  insignificant 
both  in  the  axial  and  the  radial  directions  mainly  because  of  the  high 
annulus  size  of  the  cathode  flow  channel.  Fig.  9  shows  that  the  pres¬ 
sure  drop  in  the  axial  direction,  however  insignificant,  increases  as 
the  cell  voltage  is  increased.  This  happens  as  the  volumetric  flow 
rate  of  air  decreases  down  the  channel  because  of  higher  consump¬ 
tion  of  02  at  a  lower  cell  voltage.  The  velocity  in  the  cathode  channel 
is  about  50  times  lower  than  that  of  the  anode  channel  and  also  the 
molar  flux  through  the  cathode  channel  keeps  decreasing  down  the 
channel  unlike  the  anode  channel  where  the  molar  flux  remains  the 
same  at  steady  state. 

Fig.  10  shows  the  concentration  profile  of  O2  in  the  cathode.  In 
Fig.  10,  ‘radial  distance  =  0’  is  the  cathode-electrolyte  interface  and 
‘radial  distanced’  is  the  cathode-flow  channel  interface.  Because 
of  a  very  thin  configuration  of  the  cathode,  negligible  drop  in  the 
concentration  is  observed  in  the  radial  direction  even  at  lower  volt¬ 
ages.  At  0.97  V  and  0.85  V,  the  concentration  is  found  to  be  healthy 
at  the  TPB  interface  all  along  the  cell;  however,  the  concentration 
approaches  zero  at  the  TPB  towards  the  air  exit  of  the  cell  when 
the  cell  terminal  voltage  is  0.65  V.  This  leads  the  cell  to  a  limiting 
current  density. 
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Fig.  9.  Pressure  profile  in  axial  direction  in  cathode  channel. 


4.  Effects  of  changes  in  operating  conditions 


4  A.  Flow  of  H2 


Fig.  11  shows  the  effects  of  changes  in  H2  flowrate  as  the  cell  tem¬ 
perature  is  changed.  At  a  given  temperature,  as  the  flowrate  of  H2 
is  increased,  the  current  produced  at  a  particular  voltage  increases. 
Flowever,  the  increase  in  the  cell  power  with  increase  in  FI2  flowrate 
becomes  more  pronounced  as  the  temperature  increases  more.  A 
maximum  of  10%  increase  in  current  is  observed  at  700  °C;  however, 
the  corresponding  maximum  is  about  60%  at  850  °C.  At  850  °C,  as 
the  flowrate  of  H2  is  changed  from  31  ml  min-1  to  51  ml  min-1 ,  the 
cell  comes  out  of  the  concentration  limited  region.  This  observation 
suggests  that  if  the  cell  temperature  is  increased,  it  is  beneficial  to 
increase  the  flowrate  of  FI2.  However,  the  utilization  factor  of  H2 
may  not  necessarily  increase. 


4.2.  Operating  pressure 

An  increase  in  the  cell  operating  pressure  typically  improves  the 
current,  but  such  increase  may  strongly  depend  on  the  concentra- 


Fig.  8.  Profile  of  02  concentration  in  cathode  channel. 


Fig.  10.  Profile  of  02  concentration  in  cathode. 
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Current(A) 

Fig.  11.  Effects  of  change  in  H2  flowrate  on  the  cell  I-V  characteristics. 

tion  polarization  of  the  cell.  Fig.  12  shows  the  effect  of  changes 
in  the  inlet  pressure  of  the  cell  on  its  performance  at  850  °C  and 
35  ml  min-1  H2  flow.  The  inlet  flowrates  of  both  the  channels  are 
kept  unchanged,  but  the  inlet  pressure  of  both  the  channels  are 
changed  simultaneously.  The  base  case  inlet  pressure  of  the  cell  is 
slightly  above  atmospheric  pressure.  At  this  inlet  pressure,  the  cell 
reaches  the  concentration  limited  region  as  the  voltage  falls  below 
0.6  V.  When  the  inlet  pressure  is  increased  to  1.2  atm,  the  cell  is 
no  more  concentration  limited  in  the  operating  region.  As  the  cell 
pressure  is  increased  further,  not  much  improvement  is  observed  in 
the  performance.  Therefore  the  inlet  pressure  is  not  increased  any 
further.  This  observation  is  in  agreement  with  that  of  Hussain  et  al. 
[1  ].  The  actual  improvements  in  the  performance  of  a  given  cell  due 
to  an  increase  in  the  pressure  truly  depends  on  the  operating  con¬ 
ditions  of  the  cell  and  its  hardware  design  including  the  materials 
used.  If  the  cell  has  to  be  operated  at  a  higher  pressure,  it  needs 
to  be  designed  accordingly  resulting  in  cost  increase.  Steady-state 
simulation  can  be  used  for  a  cost-benefit  analysis  to  identify  the 
optimum  operating  pressure  of  the  cell. 


Current  (A) 


Fig.  13.  Change  in  the  I-V  characteristics  of  the  cell  due  to  change  in  the  cell  oper¬ 
ating  temperature. 


4.3.  Cell  temperature 

Fig.  13  shows  the  change  in  the  current  produced  by  the  cell 
due  to  changes  in  the  temperature.  H2  flowrate  is  kept  high 
(51  ml  min-1)  for  this  study  so  that  the  cell  is  not  concentration 
limited  in  the  entire  region.  As  the  temperature  is  increased  from 
700  °C  to  750  °C,  the  improvement  in  the  performance  is  significant. 
The  current  at  0.6  V  increases  by  40%.  However  as  the  tempera¬ 
ture  is  increased  to  800  °C,  the  improvement  becomes  lower.  An 
improvement  of  about  20%  is  observed  at  0.6  V.  As  the  tempera¬ 
ture  is  increased  to  850  °C,  improvement  is  only  observed  at  lower 
cell  potential.  At  0.6  V,  the  improvement  is  only  about  5%.  Until  a 
cell  voltage  of  0.85  V,  it  is  found  that  the  cell  produces  higher  cur¬ 
rent  at  800  °C  compared  to  850  °C.  The  loss  in  the  performance  due 
to  an  increase  in  the  temperature  beyond  a  certain  value,  espe¬ 
cially  close  to  OCP,  happens  because  of  a  number  of  interacting  and 
opposing  phenomena  taking  place  in  the  cell  simultaneously.  For 
a  SOFC  with  YSZ  as  its  electrolyte,  significant  contribution  towards 
the  ohmic  loss  comes  from  the  electrolyte.  If  the  cell  temperature  is 
increased,  the  conductivity  of  the  electrolyte  improves.  It  eventually 
flattens  out  as  the  temperature  is  increased  further.  The  activa¬ 
tion  overpotential  of  both  the  electrodes  decrease  with  increase 
in  temperature.  The  diffusivity  of  the  gases  also  increases  due  to 
an  increase  in  the  temperature.  However,  the  ohmic  losses  gets 
adversely  affected  with  increase  in  the  temperature.  Also,  from  the 
ideal  gas  law,  it  can  be  observed  that  as  the  inlet  temperature  is 
increased  for  an  isothermal  cell  at  a  given  inlet  pressure,  the  con¬ 
centrations  of  the  reactants  decrease.  Therefore  the  improvement 
in  the  cell  performance  due  to  an  increase  in  the  temperature  is 
not  monotonic,  especially,  close  to  OCP.  The  effects  of  the  change 
in  the  operating  temperatures  may  vary  based  on  the  materials  of 
construction,  cell  dimensions,  and  the  control  philosophy  of  the 
cell  that  may  affect  the  spatial  variation  of  temperature  in  the  cell. 
The  results  presented  here  are  valid  for  the  modeled  industrial 
cell  which  is  of  thin-electrolyte  configuration  and  where  the  cell  is 
kept  inside  a  furnace,  its  temperature  being  maintained  by  a  feed¬ 
back  controller.  A  higher  operating  temperature  than  the  base  case 
means  higher  parasitic  losses.  In  addition,  if  the  improvement  in 
the  performance  due  to  an  increase  in  the  temperature  is  not  sub¬ 
stantially  high,  then  the  overall  cell  efficiency  may  in  fact  decrease 
with  the  increase  in  the  temperature  beyond  some  value.  A  systems 


Fig.  12.  Effects  of  cell  pressure  on  the  cell  polarization  curve. 
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level  optimization  study  may  help  to  identify  such  an  optimum 
temperature  for  a  given  cell. 

5.  Effects  of  changes  in  cell  characteristic  parameters 

5.1  Porosity 

Porosity  of  the  electrodes  play  an  important  role  in  determin¬ 
ing  the  cell  performance  along  with  maintaining  the  structural 
integrity.  An  increase  in  the  porosity  increases  the  effective  dif- 
fusivity  of  the  gaseous  species.  The  quantity  of  solid  decreases 
with  increase  in  the  porosity,  resulting  in  increased  ohmic  resis¬ 
tance.  These  opposing  effects  are  reflected  in  the  cell  performance 
shown  in  Fig.  14.  The  figure  shows  the  effect  of  changes  in  anode 
porosity  at  850  °C  and  35  ml  min-1  H2  flow.  As  the  porosity  of  the 
anode  is  decreased  below  the  base  case  porosity,  the  current  pro¬ 
duced  comes  down.  However,  as  the  porosity  is  increased  to  0.4,  an 
increase  in  the  cell  performance,  although  insignificant,  is  observed. 
As  the  porosity  is  increased  further  to  0.45,  it  results  in  a  loss  in  the 
performance.  The  inset  of  Fig.  14  makes  the  differences  very  clear. 
It  can  be  seen  that  the  optimum  porosity  of  the  existing  cell  for  this 
operating  condition  is  between  0.4  and  0.45. 

A  different  picture  is  obtained  when  the  porosity  of  the  cathode 
is  changed.  Fig.  15  shows  that  the  performance  decreases  as  the 
porosity  is  increased  by  0.1  from  the  base  case  porosity.  When  the 
cathode  porosity  is  decreased  to  0.25,  the  cell  performance  actu¬ 
ally  improves.  But  with  further  decrease  in  the  porosity,  a  loss  in 
the  performance  is  observed.  The  inset  of  Fig.  15  highlights  the  fall 
in  the  performance  with  further  decrease  in  the  cathode  porosity. 
The  effects  in  the  cathode  due  to  the  changes  in  the  porosity  are 
quite  different  from  that  of  the  anode.  The  main  reason  is  that  the 
cell  is  anode  supported  and  the  cathode  thickness  is  significantly 
smaller  than  that  of  the  anode.  In  our  previous  work  [19],  it  can  be 
observed  that  the  current  path  length  is  longer  in  the  cathode  com¬ 
pared  to  that  in  the  anode.  Also,  the  area  available  for  flow  of  the 
electrons  is  less  in  the  cathode  because  of  the  thin  configuration 
of  the  cathode.  When  the  anode  porosity  is  decreased,  the  concen¬ 
tration  loss  grew  higher  than  the  improvement  in  the  ohmic  loss. 
In  fact,  the  ohmic  loss  in  the  anode  is  negligible  at  the  base  condi¬ 
tion.  However,  in  the  cathode,  a  decrease  in  the  porosity  resulted 
in  an  decrease  of  the  ohmic  loss.  This  decrease  in  the  ohmic  loss 
is  higher  than  the  increase  in  the  concentration  loss.  For  the  given 


cell. 

operating  conditions,  the  optimum  porosity  of  cathode  is  observed 
to  be  between  0.2  and  0.25. 

5.2.  Effective  diffusivity 

A  study  is  performed  to  observe  the  effects  of  the  effective  dif¬ 
fusivity  on  the  cell.  No  specific  parameters  such  as  porosity  and 
tortuousity  are  changed  for  this  study.  Instead,  the  base  case  value 
of  the  effective  diffusivity  is  changed  to  see  its  effects  on  the  cell 
performance.  Fig.  16  shows  the  effects  of  change  in  the  effective  dif¬ 
fusivity  of  both  H2  and  H20.  As  the  effective  diffusivity  is  decreased 
to  50%  of  the  base  case  value,  a  loss  in  the  performance  in  the  entire 
region  is  observed.  However,  when  it  is  increased  by  50%  than  the 
base  case  value,  no  significant  gain  in  the  performance  could  be 
achieved.  The  reason  is  comparably  insignificant  concentration  loss 
in  the  anode  in  the  base  case  condition.  The  study  shows  that  con¬ 
siderable  improvement  in  the  performance  is  not  possible  just  by 
changing  the  effective  diffusivity  or  by  improving  the  concentration 
loss  in  the  anode. 


cell. 
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cell. 

In  the  base  case  condition,  the  cell  leads  to  the  limiting  current 
density  as  the  cell  voltage  approaches  0.55  V.  The  main  reason  is 
the  lower  effective  diffusivity  of  O2  in  the  cathode.  As  the  effec¬ 
tive  diffusivity  of  02  is  increased  by  50%  of  the  base  case  value,  the 
cell  comes  out  of  the  concentration  limited  region.  This  improve¬ 
ment  is  shown  in  Fig.  17.  But  as  the  effective  diffusivity  is  reduced 
to  50%  of  the  base  case  value,  a  significant  reduction  in  the  cell 
performance  is  observed.  The  cell  reaches  the  concentration  lim¬ 
ited  region  at  a  cell  voltage  of  about  0.75  V,  which  is  much  higher 
than  the  corresponding  base  case  voltage.  The  effects  of  the  reduc¬ 
tion  in  the  effective  diffusivity  is  significant  in  the  cathode  because 
the  base  case  value  of  the  effective  diffusivity  of  02  in  the  cath¬ 
ode  is  much  lower  than  that  of  H2  and  H20  in  the  anode.  This 
study  suggests  that  efforts  should  be  made  to  lower  the  concentra¬ 
tion  overpotential  in  the  cathode  in  order  to  avoid  limiting  current 
density  especially  when  higher  currents  are  produced  from  the 
cell. 

5.3.  Reaction  rate 

A  study  is  done  to  observe  the  effects  of  the  reaction  rate  on  the 
overall  performance  of  the  cell.  Since  the  cell  is  believed  to  follow 
the  Butler-Volmer  kinetics,  reaction  rates  are  modified  by  chang¬ 
ing  the  pre-exponential  factors  in  the  equation  for  the  exchange 
current  density  [19].  When  the  anodic  reaction  rate  is  modified 
by  ±20%,  the  performance  of  the  cell  remains  almost  unaffected. 
This  shows  that  the  anodic  activation  overpotential  is  negligible 
for  this  cell.  Therefore  these  results  are  not  shown.  Fig.  18  shows 
the  effects  of  the  change  in  the  cathodic  reaction  rate.  When  the 
value  of  kca  is  reduced  by  10%  compared  to  the  base  case,  the  cell 
performance  does  not  change  appreciably  up  to  a  voltage  of  0.9 
V.  The  difference  in  performance  becomes  wider  as  the  cell  volt¬ 
age  is  further  reduced.  At  0.6  V,  the  current  produced  is  lower  by 
6.25%.  Similar  results  are  seen  as  kca  value  is  increased  by  10%.  How- 
ever,  the  improvement  in  the  performance  is  lower  compared  to  the 
reduction  in  the  performance  for  a  corresponding  decrease  in  the 
pre-exponential  factor.  The  reason  is  that  higher  losses  (concen¬ 
tration  and  ohmic  losses)  result  as  the  cell  produces  more  current 
due  to  an  increase  in  the  pre-exponential  factor  in  the  expression 
for  cathodic  exchange  current  density.  At  0.6  V,  the  current  pro¬ 
duced  increases  by  5.9%.  Beyond  0.6  V,  the  percentage  change  in 
the  cell  performance  compared  to  the  base  case  becomes  almost 
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Fig.  18.  Effects  of  the  change  in  cathode  reaction  rate  on  the  I-V  characteristics  of 
the  cell. 

constant  for  both  the  cases  as  the  concentration  loss  starts  domi¬ 
nating. 

6.  Effects  of  changes  in  cell  design  parameters 

6.1.  Thickness  of  the  anode 

In  the  following  section,  a  study  is  presented  by  decreasing  the 
thickness  of  the  electrodes  and  the  electrolyte.  If  the  thickness  of 
the  anode  is  decreased  for  an  anode-supported  tubular  cell,  the 
concentration  overpotential  at  the  anode  decreases  because  of  a 
lower  diffusion  path  length.  A  smaller  thickness  of  the  anode  results 
in  a  smaller  inner  radius  of  the  electrolyte  and  the  cathode.  The 
smaller  inner  radius  of  the  cathode  results  in  a  shorter  current  path 
length  because  of  the  location  of  the  current  collector.  The  ohmic 
loss  in  the  anode  also  reduces.  However,  a  smaller  thickness  of  the 
anode  will  result  in  a  smaller  active  area  for  the  electrochemical 
reactions  at  the  TPB  of  both  the  electrodes  increasing  the  overall 
activation  overpotential.  It  can  be  seen  in  Fig.  19  that  as  the  anode 
thickness  is  decreased,  there  is  an  improvement  in  the  overall  cell 
performance  compared  to  the  base  case  till  a  voltage  of  0.6  V.  But 
the  limiting  current  density  of  the  cell  is  lower  than  the  base  case.  As 
the  circumferential  area  available  for  the  flux  of  the  gaseous  species 
to  the  TPB  decreases  with  a  decrease  in  the  anode  thickness,  the 
cell  goes  to  a  limiting  current  region.  When  the  anode  thickness 
is  increased  by  30%  of  the  base  case  value,  the  current  produced 
from  the  cell  decreases.  However,  because  of  the  reasons  mentioned 
before,  the  cell  is  not  concentration  limited  till  a  voltage  of  0.5  V 
compared  to  0.55  V  in  the  base  case.  In  summary,  it  can  be  said  that 
the  decrease  in  the  anode  thickness  improves  the  cell  performance 
for  most  of  the  polarization  curve,  but  the  limiting  current  density 
gets  decreased.  The  reverse  is  observed  as  the  thickness  of  the  anode 
is  increased. 

6.2.  Thickness  of  the  electrolyte 

As  the  thickness  of  the  electrolyte  is  decreased,  the  ohmic  loss 
that  occurs  in  the  electrolyte  is  reduced.  Also  the  current  path 
length  in  the  cathode  decreases.  There  is  no  change  in  the  active 
area  for  the  anodic  reaction,  but  the  active  area  for  the  cathodic 
reaction  gets  reduced.  In  Fig.  20,  it  is  seen  that  as  the  thickness  of 
the  electrolyte  is  reduced  by  30%,  the  cell  performance  improves. 
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Fig.  19.  I-V  characteristics  due  to  change  in  anode  thickness. 


However,  as  the  thickness  of  the  electrolyte  is  already  very  small 
in  the  base  case,  the  reduction  does  not  result  in  much  decrease  in 
the  inner  radius  of  the  cathode.  The  result  is  a  higher  current  den¬ 
sity  till  0.5  V.  As  the  electrolyte  thickness  is  increased  by  30%,  the 
dominating  ohmic  loss  in  the  electrolyte  decreases  the  cell  perfor¬ 
mance. 

6.3.  Thickness  of  the  cathode 

Fig.  21  shows  that  change  in  the  thickness  of  the  cathode  has 
a  very  strong  effect  on  the  cell  performance.  When  the  cathode 
thickness  is  decreased,  the  area  for  flow  of  the  electrons  from  the 
cathode  current  collector  to  the  TPB  gets  decreased.  This  adds  on 
to  the  existing  significant  ohmic  loss  in  the  cathode.  Even  though 
the  diffusion  path  length  decreases  due  to  the  decrease  in  cathode 
thickness,  the  overall  effect  is  a  loss  in  the  performance.  The  cur¬ 
rent  from  the  cell  keeps  decreasing  as  the  cell  voltage  is  lowered. 
With  the  increase  in  the  thickness,  the  current  produced  increases 
because  of  a  lower  ohmic  loss  in  the  cathode.  However  because  of 
the  increased  diffusion  length,  the  concentration  limitations  start 
showing  up  at  a  higher  voltage  than  the  base  case. 


7.  Dimensional  optimization  study 

In  order  to  identify  whether  an  optimum  combination  for  the 
dimensions  of  the  anode  channel  radius,  cell  length,  and  the  cath¬ 
ode  channel  annulus  exists,  a  dimensional  study  is  done.  In  this 
study,  the  thickness  of  the  electrodes  and  the  electrolyte  are  kept 
constant  at  the  base  case  value.  The  radius  of  the  anode  channel 
is  varied.  As  the  volume  of  the  solid  is  kept  constant,  a  new  set 
of  values  for  the  cell  length  are  generated.  As  the  total  volume  of 
the  cell  (excluding  the  insulator)  is  also  kept  constant,  the  annu¬ 
lus  size  is  changed  depending  on  the  new  cell  length  and  the  new 
anode  channel  radius.  The  lowest  radius  considered  in  this  study 
is  one  that  can  be  manufactured  with  the  commercially  available 
processes.  It  is  observed  in  Fig.  22  that  there  is  an  optimum  power 
density  of  the  cell  which  does  not  correspond  to  the  minimum 
value  of  the  anode  channel  radius.  It  is  also  observed  that  for  any 
radius  larger  than  the  optimum,  the  power  density  monotonically 
decreases. 


Fig.  20.  I-V  characteristics  due  to  change  in  electrolyte  thickness. 


Fig.  22.  Power  density  vs.  length  and  anode  channel  radius. 
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Fig.  23.  Pressure  drop  and  ohmic  polarization  with  change  in  anode  channel  radius. 

The  competing  loss  mechanisms  in  the  cell  are  depicted  in 
Fig.  23.  It  is  observed  that  the  pressure  drop  across  both  the  anode 
and  cathode  channels  increase  with  a  decrease  of  the  anode  channel 
radius.  This  causes  an  overall  rise  in  the  concentration  loss  but  the 
ohmic  loss  decreases  with  a  decrease  in  the  anode  channel  radius. 
The  steeper  rise  in  the  pressure  drop  in  the  anode  channel  can  be 
attributed  to  the  fact  that  the  velocity  in  the  anode  channel,  even 
at  the  base  case,  is  very  high.  Further  decrease  in  the  radius  causes 
this  velocity  to  rise  higher.  Furthermore,  more  water  is  generated 
with  the  increase  in  the  current.  The  higher  concentration  of  water 
vapor  causes  the  pressure  drop  to  increase  because  of  the  higher 
viscosity  of  water  vapor  than  that  of  FI2.  The  activation  overpoten¬ 
tial  also  gets  affected  because  of  the  change  in  the  anode  channel 
radius.  Because  of  these  competing  loss  mechanisms,  an  optimum 
power  density  can  be  achieved  at  an  inner  channel  radius  that  can 
be  manufactured  commercially.  Flowever,  a  nonlinear  constrained 
optimization  needs  to  be  done  to  reach  the  optimum  set  of  values 
considering  more  objective  functions  in  a  wide  operating  range. 
Such  a  study  is  presented  elsewhere  [19]. 

8.  Conclusions 

The  study  in  this  paper  shows  that  the  radial  change  in  the 
species  concentration  inside  the  anode  flow  channel  can  become 
significant  when  the  cell  produces  higher  current.  Similar  conclu¬ 
sions  are  drawn  for  the  H2  and  H20  concentrations  in  the  anode. 
Even  though  the  pressure  drop  in  the  anode  flow  channel  is  insignif¬ 
icant,  it  may  become  significant  when  the  cell  produces  more  water. 
As  the  pressure  in  the  cathode  channel  is  near  atmospheric,  02 
concentration  is  found  to  approach  zero  at  the  channel-cathode 
interface  at  high  current  densities.  These  variations,  therefore,  need 
to  be  considered  for  designing  cells  producing  higher  currents  and 
for  optimization  studies.  Pressure  drop  in  the  cathode  channel  is 
found  to  be  insignificant  and  lesser  than  that  in  the  anode  chan¬ 
nel.  In  addition,  the  pressure  drop  is  less  at  high  current  densities 
as  02  gets  consumed.  Because  of  the  very  thin  configuration  of 
the  cathode,  the  concentration  loss  in  the  radial  direction  inside 
the  cathode  is  not  significant  except  at  higher  current  densities. 
An  increase  in  the  flowrate  of  FI2  boosts  the  cell  performance  in 
the  operating  range  considered  in  this  study.  The  improvement  is 
higher  when  the  temperature  is  higher.  Flowever  the  utilization 


factor  of  H2  can  be  adversely  affected.  Beyond  a  certain  pressure, 
no  significant  improvement  in  the  performance  is  observed.  The 
increase  in  the  performance  due  to  an  increase  in  the  temperature 
is  not  monotonic,  especially  close  to  OCP.  The  best  operating  condi¬ 
tions  for  the  cell  can  be  identified  by  an  optimization  study  at  the 
systems  level  by  considering  the  net  power  output  from  the  cell  and 
a  cost-benefit  analysis.  The  optimum  porosity  for  the  given  operat- 
>  ing  condition  is  found  to  be  0.4-0.45  for  the  anode  and  0.2-0.25  for 
^  the  cathode.  It  is  also  observed  that  the  cell  performance  cannot  be 

£  improved  by  decreasing  the  concentration  loss  in  the  anode.  FIow- 

■JJ  ever,  a  significant  improvement  in  the  performance  is  possible  by 
c  decreasing  the  concentration  loss  in  the  cathode.  This  is  important, 
O  in  particular,  when  the  current  density  is  high.  The  anodic  activa¬ 
tion  loss  is  found  to  be  insignificant.  Flowever,  the  performance  can 
be  strongly  affected  by  changing  the  activation  loss  in  the  cath¬ 
ode  till  the  concentration  loss  starts  dominating.  A  decrease  in  the 
anode  thickness  for  this  anode-supported  cell  results  in  an  overall 
improvement  in  the  performance  but  the  limiting  current  density 
is  decreased.  A  decrease  in  the  electrolyte  thickness  boosts  the  cell 
performance  mainly  because  of  a  decrease  in  the  ohmic  loss.  A  very 
strong  effect  of  the  thickness  of  the  cathode  is  observed.  When  the 
cathode  thickness  is  decreased,  there  is  an  overall  loss  in  the  per¬ 
formance.  An  increase  in  the  cathode  thickness  reduces  the  ohmic 
loss,  but  leads  the  cell  quickly  to  a  concentration  limited  region.  The 
dimensional  sensitivity  study  shows  that  it  is  possible  to  obtain  a 
set  of  values  for  the  anode  channel  radius,  the  length  of  the  cell,  and 
the  annulus  size  that  can  result  in  an  optimum  power  density.  The 
anode  channel  radius  that  results  in  the  optimum  power  density 
can  be  manufactured  commercially. 
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